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Abstract. Experimental results for congested pedestrian traffic are pre- 
sented. For data analysis we apply a method providing measurements on 
an individual scale. The resulting velocity-density relation shows a co- 
existence of moving and stopping states revealing the complex structure 
of pedestrian fundamental diagrams and supporting new insights into 
the characteristics of pedestrian congestions. Furthermore we introduce 
a model similar to event driven approaches. The velocity-density relation 
as well as the phase separation is reproduced. Variation of the parame- 
ter distribution indicates that the diversity of pedestrians is crucial for 
phase separation. 
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1 Introduction 

The velocity-density relation is one of the most important characteristics for 
the transport properties of any traffic system. For pedestrian traffic there is 
currently no consensus even about the principle shape of this relation which is 
reflected e.g. in conflicting recommendations in various handbooks and guidelines 
[1]. Discrepancies occur in particular in the high-density regime vifhich is also 
the most relevant for applications in safety analysis like evacuations or mass 
events. At high densities stop-and-go waves occur indicating overcrowding and 
potentially initiating dangerous situations due to stumbling etc. However, the 
densities where the flow breaks down due to congestion ranges from densities 
of /Cmax = 3.8 m~^ to /Omax = 10 m~^ [T] . This large variation in values for 
Pmax reported in the literature is partly due to insufficient methods of data 
capturing and data analysis. In previous experimental studies, different kinds of 
measurement methods are used and often a mixture of time and space averages 
are realized. Especially in the case of spatial and temporal inhomogeneities the 
choice of the measurement method and the type of averaging have a substantial 
influence on the results [2]. 

Up to now, congested states in pedestrian dynamics have not been analyzed 
in much detail. This is in contrast to vehicular traffic where the congested phase 
is well-investigated, both empirically and theoretically 
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In this contribution we show that even improved classieal measurement meth- 
ods using high precision trajectories but basing on mean values of density and 
velocity fail to resolve important characteristics of congested states. For a thor- 
ough analysis of pedestrian congestion we apply a new method enabling mea- 
surements on the scale of single pedestrians. 



2 Experimental data 

For our investigation we use data from experiments performed in 2006 in the 
wardroom of Bergische Kaserne Diisseldorf with a test group of up to iV ~ 70 
soldiers. The length of the circular system was about 26 m, with a Z = 4 m long 
measurement section. Detailed information about the experimental setup and 
data capturing providing trajectories of high accuracy (|xorr| < 0.02 m) is given 
in PI5] . 
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Fig. 1. Trajectories for the runs with N = 39, 56, 62 and 70 (left to right, top to 
bottom). With increasing density the occurrence of stop-and-go waves accumu- 
late. 
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In Fig. [1] the x-componcnt of trajectories is plotted against time. For the 
extraction of the trajectories, the pedestrians' heads were marked and tracked. 
Backward movement leading to negative velocities is caused by head movement 
of the pedestrians during a standstill. Inhomogeneities in the trajectories increase 
with increasing density. As in vehicular traffic, jam waves propagating opposite 
to the movement direction (upstream) occur at higher densities. 

Stopping is first observed during the runs with A'^ = 45 pedestrians, at 70 
pedestrians they can hardly move forward. Macroscopically one observes sepa- 
ration into a stopping area and an area where pedestrians walk slowly. In the 
following we analyze how macroscopic measurements blur this phase separation 
and apply a technique introduced in [B] enabling a measurement of the funda- 
mental diagram on a 'microscopic' scale. 

2.1 Macroscopic measurement 

Speed Vi of pedestrian i and the associated density pi are calculated using the 
entrance and exit times t™ and into and out of a measurement section of 
length /„! = 2 m, 

= = ^J-^. r Pit) dt with p{t) = . (1) 




Fig. 2. Velocity-density relation using an improved macroscopic measurement 
method. There is no indication of phase separation since stopping states {v ~ 0) 
occur only at large densities where no moving states {v > 0) are observable. 



The speed Vi is a mean value over the space-time interval Ati = — 
and Ax = Im- By integration over the instantaneous density p{t) the density is 
assigned to the same space-time interval. To reduce the fluctuations of p{t) we 
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use the quantity Oi{t), introduced in [7], which measures the fraction of space 
between pedestrians i and i + I that is inside the measurement area. 

Results of the macroscopic measurement method are shown in Fig. [2l In 
comparison to method B introduced in [J] (see Fig. 6 of [2] which uses data 
based on the same trajectories as this study) the scatter of the data is reduced 
due to an improved density definition and a better assignment of density and 
velocity. However the resulting velocity-density relation does not allow to identify 
phase-separated states although these are clearly visible in the trajectories. 

2.2 Microscopic measurement 

To identify phase separated states we determine the velocity-density relation 
on the scale of single pedestrians. This can be achieved by the Voronoi density 
method [B]. In one dimension a Voronoi cell is bounded by the midpoints Zi = 
(xi+i +Xi)/2 of the pedestrian positions Xi and x^+i. With the length Li = 
Zi — Zi-i of the Voronoi cell corresponding to pedestrian i and At = 0.5 s we 
define the instantaneous velocity and density by 

v'M = ^^(^ + -^V2) -..ft-^V2) ^ r l/L : X 6 [z. 

^ At ' ^ 1 : otherwise ^ ^ 




Fig. 3. Left: Velocity-density relation on an individual scale. At each density 
p' > 1.5 m~^ both stopping and moving states are observable. Right: Probabil- 
ity distribution of the velocities for different density intervals. The double peak 
structure indicates the coexistence of moving and stopping states. The height of 
the stopping peak increases with increasing density. 



The fundamental diagram based on the Voronoi method is shown on the left 
side of Fig. [3] Regular stops occur at densities higher than 1.5 m~^. On the right 
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side of Fig. |3]thc distribution of the velocities for fixed densities from 1.8 m"-'^ 
to 2.6 ni~^ are shown. There is a continuous change from a single peak near 
V = 0.15 m/s, to two peaks, to a single peak near v = m/s. The right peak 
represents the moving phase, whereas the left peak represents the stopping phase. 
At densities around 2.2 m^^ these peaks coexist, indicating phase separation into 
a flowing and a jammed phase. 

2.3 Phase separation in vehicular traffic 

In highway traffic, phase separation into moving and stopping phases typically 
occurs when the outflow from a jam is reduced compared to maximal possible 
flow in the system. Related phenomena are hysteresis and a non-unique fun- 
damental diagram. At intermediate densities two different flow values can be 
realized. The larger flow, corresponding to a homogeneous state, is metastable 
and breaks down due to fluctuations or perturbations (capacity drop) . The origin 
of the reduced jam outflow is usually ascribed to the so-called slow-to-start be- 
haviour (see |3l4j and references therein), i.e. an delayed acceleration of stopped 
vehicles due to the loss of attention of the drivers etc. 

The structure of the phase-separated states in vehicular traffic is different 
from the ones observed here. For vehicle traffic the stopping phase corresponds 
to a jam of maximal density whereas in the moving phase the flow corresponds to 
the maximal stable flow, i.e. all vehicles in the moving phase move at their desired 
speed. This scenario is density-independent as increasing the global density will 
only increase the length of the stopping region without reducing the average 
velocity in the free flow regime. The probability distribution of the velocities (in 
a periodic system) shows a similar behaviour to that observed in Fig. |31 The 
position of the free flow peak in the case of vehicular traffic is independent of 
the density. 

The behaviour observed here for pedestrian dynamics differs slightly from 
that described above. The main difference concerns the properties of the moving 
regime. Here the observed average velocities are much smaller than the free 
walking speeds. Therefore the two regimes observed in the phase separated state 
are better characterized as "stopping" and "slow moving" regimes. 

Further empirical studies are necessary to clarify the origin of these differ- 
ences. One possible reason are the different acceleration properties of vehicles 
and pedestrians as well as anticipation effects. It also remains to be seen whether 
in pedestrian systems phenomena like hysteresis can be observed. 

3 Modeling 

3.1 Adaptive Velocity Model 

In this section we introduce the adaptive velocity model, which is based on an 
event driven approach [Hj . A pedestrian can be in different states which determine 
the velocity. A change between these states is called event. The model was derived 
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from force-based models, where the dynamics of pedestrians are given by the 
following system of coupled differential equations 

m,^=i^, with i^, = J^drv ^ ^rcp ^^^^ ^ = t'. , (3) 

where Fi is the force acting on pedestrian i. The mass is denoted by m^, the 
velocity by Vi and the current position by Xi. Fi is split into a repulsive force 
F™'' and a driving force Ff"^. The dynamics is regulated by the interrelation 
between driving and repulsive forces. In our approach the role of repulsive forces 
are replaced by events. The driving force is defined as 

n 

where v'^ is the desired speed of a pedestrian and r the relaxation time of the 
velocity. By solving the differential equation 

at \ n J 

the velocity function is obtained. This is shown in Fig. |4] together with the 
parameters governing the pedestrians' movement. In this model pedestrians are 




Fig. 4. Left: Connection between the required space d, the step length st and 
the safety distance /3. Right: The adaptive velocity: acceleration until idcci than 
deceleration until iacci, again acceleration until tdcc2 and so on. 



treated as bodies with diameter d [S]. The diameter depends linearly on the 
current velocity and is equal to the step length st in addition to the safety 
distance /3 

=st,(<)+A(t) (6) 

Step length and safety distance are introduced to define the rules for the 
dynamics of the system. We determine the model parameters from empirical 
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data which allows to judge the adequacy of the rules. Based on [TU] the step 
length is a linear function of the current velocity with following parameters: 

sti{t) = 0.235m + 0.302[s] v^{t). (7) 

The required quantities for the safety distance can be specified through em- 
pirical data of the fundamental diagram di = 1/p = 0.36 + 1.06 v, see [7]. With 
these experimental results the previous equations can be summarized to 

ft (t) =d,{t)~ sU {t)^a, + h [t) (8) 

with ai = 0.125 m and hi = 0.758 s. No free model parameter remain with 
these specifications. In the following we describe the rules for the movement. A 
pedestrian accelerates to the desired velocity vf until the distance Axi^i+i to the 
pedestrian in front is smaller than the safety distance. From this time on, he/she 
decelerates until the distance is larger than the safety distance. To guarantee a 
minimal volume exclusion, case "collision" is included, in which the pedestrians 
are too close to each other and have to stop. Via Axi^i+i, di and ft; the velocity 
function for the states deceleration (dec), acceleration (acc.) and collision (coll.) 
can be defined, see Eq. O 

Wdoc exp (-■^^^^) for Axis+i - Si{t) < (dec.) 
v,{t) = { v^{l - exp (-^) for At,- ,+i - S^{t) > (acc.) (9) 

for Ax.^+i^d,{t) <~/3,{t)/2 (coll.) 

where Si{t) ~ {di{t) + di^i{t))/2 is the distance between the centers of both 
pedestrians. The current velocity Vi(t) of an pedestrian i depends on his/her 
state, ^dec denominates the point in time where a change from acceleration to 
deceleration takes place. Conversely tacc is the change from deceleration to ac- 
celeration. Vdcc = v{tdec) and Vacc = i'(^acc) are defined accordingly. At the 
beginning to = s with a change from acceleration to deceleration a new calcu- 
lation of to is necessary: 

to = tacc + In fl - ^ exp f ih^l^lf^) ) (10) 



The discreteness of the time step could lead to configurations where overlap- 
ping occurs. To ensure good computational performance for high densities, no 
events are explicitly calculated. Instead in each time step, it is checked whether 
an event has taken place and tdec, ^acc or tcoU are set to t accordingly. To avoid 
too large interpenetration of pedestrians and to implement a reaction time in 
a realistic size we choose At = 0.05 s. To guarantee a parallel update a recur- 
sive procedure is necessary: Each person is advanced one time step according to 
Eq. [9] If after this step a pedestrian is in a different state because of the new 
distance to the pedestrian in front, the velocity is set according to this state. 
Then the state of the next following person is reexamined. If the state is still 
valid the update is completed. Otherwise, the velocity is calculated again. 
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3.2 Model validation and influence of individual diff'erences 

In the following we face model results with experimental data and study how 
the distribution of individual parameter influences the phase separation. For all 
simulations the desired velocity is normal distributed w,^ ~ A/'(/x, cr^) with average 
II = 1.24 m/s and variance cr^ — 0.05. Fig. [5] and Fig. [5] show the simulation 
results for two different choices of parameter distributions. 
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Fig. 5. Validation of the modeled fundamental diagram and trajectories with 
same parameter Ui, hi and for all pedestrians Left: Comparison of fundamental 
diagrams of modeled and empirical data. Middle: Trajectories for = 62 
(model). Right: Trajectories for N = 7Q (model). 




Fig. 6. Validation of the modeled fundamental diagram and trajectories, with 
normal distributed individual parameter. Left: Comparison of fundamental dia- 
grams of modeled and empirical data. Middle: Trajectories for A^ = 62 (model). 
Right: Trajectories for A^ = 70 (model). 



The model yields the right macroscopic relation between velocity and den- 
sity even if ai,bi,sti and are the same for all pedestrians, see Fig. [S] (left). 
The trajectories display that phase separation does not appear. Even at high 
densities the movement is ordered and no stops occur. For further simulations 
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we incorporate a certain disorder by choosing the following individual parameter 
normal distributed: a, - A/'(0.125, 0.1), b, ~ A/'(0.758, 0.5) and n - A/'(1.0, 0.1). 

Variation of the personal parameters affects the scatter of the fundamental 
diagram, see Fig. [51 left. Phase separation appears in the modeled trajectories 
as in the experiment, see Fig. [5] middle and right. It is clearly visible that long 
stop phases occur by introducing distributed individual parameters. Then the 
pattern as well as the change of the pattern from iV = 62 to TV = 70 are in good 
agreement with the experimental results, see Fig. [TJ Even the phase separated 
regimes match qualitatively. However, the regimes appear more regular in the 
modeled trajectories. 

In Fig. [7] microscopic measurements of fundamental diagram and the related 
velocity distributions arc shown. Separation of phases is reproduced well. But 
the position of the peak attributed to the moving phase is not in conformance 
with the experimental data, compare Fig. [3] (right). The experimental data show 
that the peak position is independent from the density at v around 0.15 m/s. 
At the model data the position of the peak changes with increasing density. 
Measurements with different time steps show that the size of the time step 
influences the length and shape of the stop phase at high densities. But the 
density where first stops occur seems independent from the size of the time 
step. Further model analysis is necessary to study the role of the reaction time 
implemented by discrete time steps in this special type of update. Furthermore 
we will study how the change of the peak could be influenced by including a 
distribution for the step length and other variations of the distribution for the 
safety distance. 




Fig. 7. Left: Microscopic fundamental diagram of the modeled data with normal 
distributed individual parameters. Right: Distribution of v at fixed densities. 
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4 Conclusion 

Wc have investigated the congested regime of pedestrian traffic using high- 
quality empirical data based on individual trajectories. Strong evidence for phase 
separation into standing and slow moving regimes is found. The corresponding 
velocity distributions show a typical two-peak structure. The structure of the 
trajectories is well reproduced by an adaptive velocity model which is a variant 
of force-based models in continuous space. 

Future studies should clarify the origin of the differences to the phase sepa- 
rated states observed in vehicular traffic. Here phase separation into a stopping 
and a moving phase occurs such that the average velocity in the moving regime 
is independent of the total density. 
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